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In this paper we study non-standard neutrino interactions as an example of piiysics beyond tlie 
standard model using atmospheric neutrino data collected during the Super-Kamiokande 1(1996- 
2001) and 11(2003-2005) periods. We focus on flavor-changing-neutral-currents (FCNC), which 
allow neutrino flavor transitions via neutral current interactions, and effects which violate lepton 
non-universality (NU) and give rise to different neutral-current interaction-amplitudes for different 
neutrino flavors. We obtain a limit on the FCNC coupling parameter, £f^r, Ie/jtI < 1.1 x 10~^ at 
90%C.L. and various constraints on other FCNC parameters as a function of the NU coupling, See- 
We find no evidence of non-standard neutrino interactions in the Super-Kamiokande atmospheric 
data. 
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PACS numbers: 13.15.+g, 14.60.Pq, 14.60.St 



I. INTRODUCTION 

The experimental understanding of neutrino oscilla- 
tions has improved dramatically over the last ten years. 
In 1998 Super-Kamiokande ("Super-K") reported an up- 
down asymmetry in the zenith angle distribution of 
muon-like (^-like) events and concluded the distortion 
was evidence for neutrino oscillations [l|. Super-K also 
observed an oscillation signature consistent with the L / E 
(pathlength over energy) dependence predicted by to 
Vr oscillations at maximal ^23 mixing 0- Results con- 
sistent with the Vy — > v-r oscillation hypothesis have 
also been obtained by the tau appearance analysis in 
Super-K g and the long-baseline accelerator experi- 
ments K2K and MINOS @. 

However, many alternatives to neutrino oscillations 
have been proposed to explain the asymmetry of the at- 
mospheric neutrino /i-like event s amp le @ . Neutrino de- 
coherence 0,[11j neutrino decay @,[lo|, mass- varying neu- 
trinos [ll| , and CPT violation effects are among the most 
prominent. Although most of these approaches have been 
ruled out by the Super-K 0, ^ and MINOS [H data, 
non-standard neutrino interactions (NSI) with matter, 
that is interactions which are not predicted by the Stan- 
dard Model, remain viable. 

Among the many types of NSI models, we focus in 
this paper on two generic types of neutrino interac- 
tions. Flavor-changing-neutral-current (FCNC) effects 
represent neutrino interactions with fermions, /, in mat- 
ter that induce neutrino flavor change: I'a + f ^ i^p + f , 
where a and /? denote neutrino flavors. Lepton non- 
universal (NU) interactions, on the other hand, are de- 
fined by a non-universal neutral current scattering ampli- 
tude among the three flavored neutrinos. In the Standard 
Model, this amplitude is identical among the neutrinos. 
Other theories of neutrino mass, however, predict these 
kinds of interactions. For instance, NSI interactions are 
often seen in models where the neutrino mass arises from 



admixtures of isosingiet neutral heavy leptons |14| or 
from i?-parity violating supersymmetry Il5|. A summary 
of various NSI models is presented in Ref. [l6[ . Although 
NSI are predicted by various theories, the expected phe- 
nomena do not generally depend on the particular phe- 
nomenologieal model and are typically characterized by 
dependence on the neutrino energy and the surrounding 
matter density. For this reason, NSI can be explored in 
a general context using atmospheric neutrinos. While 
the current constraints on NSI come from beam-based 
experiments [sil . [3^ , atmospheric neutrinos can provide 
additional sensitivity to these interactions due to their 
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ample flux and the large amount of matter they traverse 
before detection. 

This paper discusses atmospheric neutrino oscillations 
in the context of NSI at Super-K and is organized as fol- 
lows. The data set and oscillation framework used in this 
paper are presented in See. |lT] and Sec. IIIIl respectively. 
In See. IIVI we show the results of an analysis assuming 
two-flavor i^^ Vr neutrino oscillations amidst NSI. An 
analysis using an extended three-flavor framework is per- 
formed in Sec. IVland finally, the results are summarized 
in Sec. EIIl 



II. ATMOSPHERIC NEUTRINO DATA 

Super-K is a 50 kiloton water Chercnkov detector lo- 
cated in a zinc mine in Kamioka, Japan. It is optically 
separated into an Inner Detector (ID) which is istru- 
mented with 11,146 inward facing photomultiplier tubes 
(PMTs) at full capacity, and an outer detector (OD) used 
to veto cosmic ray muons. A more detailed description 
of the detector is presented in Ref. [l7| . The run period 
of Super-K has been classified into four phases: The first 
phase corresponds to physics data taken between April 
1996 and July 2001 (SK-I). After an accident at the end 
of 2001, Super-K resumed data taking with half the num- 
ber of ID PMTs between October 2002 and October 2005 
(SK-II). The remaining two run periods are divided into 
SK-III (2006-2008), after rebuilding the ID with its full 
complement of PMTs, and SK-IV (2008-present), after 
the data acquisition system was upgraded. In this paper 
we use data from the SK-I and SK-II run periods. 



1. Classification of atmospheric neutrino data 

Atmospheric neutrino events in Super-K are di- 
vided into the following four categories: fully con- 
tained (FC), partially contained (PC), upward stopping 
muons (UPMU stopping) and upward through-going 
muons (UPMU through-going). For FC and PC events, 
event vertices are required to be within the ~22.5 kton 
fiducial volume defined by the volume inset from the ID 
walls by 2 m. An event whose particles are completely 
contained within the ID is classified as FC, while an event 
with particles exiting the ID and depositing energy into 
the outer detector (OD) is classified as PC. The PC sam- 
ple is further classified into two sub-categories, PC stop- 
ping and PC through-going. The former corresponds to 
events with a particle that stops in the OD, while in the 
case of the latter, the particle exits the OD. UPMU events 
are produced by the charged current interactions of at- 
mospheric muon neutrinos in the rock surrounding the 
detector. Muons traveling in the upward direction are 
selected to avoid contamination from cosmic ray muons. 
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The UPMU stopping sample is defined by events which 
enter from outside the detector and stop inside the ID, 
while the UPMU through-going sample is composed of 
those that enter and subsequently exit the ID. The ex- 
pected mean energy for each of the event class is, ^ 1 GeV 
for FC, - 10 GeV for PC, - 10 GeV for UPMU stopping, 
and ^ 100 GeV for UPMU through-going. Descriptions 
of the event reduction and reconstruction can be found 
in 

2. Monte Carlo simulation 

In this paper, independent 500 year Monte Carlo (MC) 
samples are used for SK-I and SK-II. The analyses use 



In this equation U is the unitary PMNS matrix [20| . 
which describes standard neutrino mixing as rotations 
among pairs of mass eigenstates parameterized by unique 
mixing angles, 9ij . Here the squared difference of the neu- 
trino masses is denoted by Am|^ , T4isw is the MSW po- 
tential in the flavor basis [2l|, Gf is the Fermi coupling 
constant, TV/ is the fermion number density in matter 
along the path of the neutrino, and the Sap represent 
the NSI coupling parameters. The non- universal cou- 
plings are represented by the flavor diagonal e and the 
FCNC interactions by the off-diagonal elements. Stan- 
dard neutrino oscillations are recovered when all of the 
£a0 = 0. Note that while the first term of the Hamil- 
tonian carries an explicit energy dependence, the second 
and third terms do not and instead are functions of the 
local matter density. For our calculations below we em- 
ploy the PREM model ^ of the Earth's density profile 
and chemical composition, where the proton to nucleon 
ratio in the mantle and core are set to be Yp = 0.497 
and 0.468, respectively [25|. In many of our calculations 
we use the average matter density along the path of the 
neutrino. 



IV. ANALYSIS WITH A TWO-FLAVOR 
HYBRID MODEL 

We consider first a model in which NSI effects in the 
i^fi — sector coexist with standard two-flavor 4-> Vt 
neutrino oscillations. In this scenario all NSI that couple 
to Ve in Eq.(IT]), e^^, are set to zero. Allowing the remain- 
ing parameters which couple to v^j_ and Vr to be non-zero 
introduces a matter-dependent effect on the oscillations 



the IIonda2006 neutrino flux [T^ and neutrino interac- 
tions are simulated using the NEUT interaction genera- 
tor [lii. 



HI. OSCILLATION FRAMEWORK 

In the following sections we will consider two separate 
NSI models derived from a more general formalism. We 
introduce first the more general framework, restricting its 
scope later to the particular NSI effects we aim to study. 
The full three-flavor Hamiltonian, Hap, governing the 
propagation of neutrinos in the presence of effects from 
NSI is 



(1) 

I 

of i^fj^ ^ Vt- Since the standard two-flavor scenario {9i2, 
013, and Am^i = 0) does not incorporate oscillations into 
i/e, there is no separate effect from the standard matter 
potential, Vmsw- Labeling this the 2- flavor hybrid models 
we can explore NSI couplings by searching for matter- 
induced distortions of standard i'^ oscillations. 



A. Formalism 

The 2-flavor hybrid model can be extracted from 
Eq.([T]) by setting Ani^i ~ 0. Following the formal- 
ism of M.C. Gonzalez- Garcia and M. Maltoni [l^], as- 
suming that neutrinos possess non-standard interactions 
with only d-quarks [H, [11], and deflning e,,T- = e and 
Err ~ £fifi = s' the survival probability in constant 
density matter is given by 

P,^_,,, = 1 - P,^^,^ = 1 - sin' 2esin2 (^^i?) • 

(2) 

The effective mixing angle, 0, and the correction factor 
to the oscillation wavelength, i?, are given by 

sin^ 26 = {sm''2e + Rlsm'2^ + 2i?osin26lsin2^) , 
R = ^l + Rl + 2i?o(cos26lcos2^ + sin26lsin2^), 

^ = ^tan-i (I) , (3) 



]^ /O \ / eee £efi EerX 

Hap = TF^Uc^j ^"^ii {U'')kp + Vmsw + v^GpNf e^t, £^,^, Sf^r 
\0 Am§J \e,r e^r ^rr) 
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where 9 is the standard two-flavor mixing angle, and ^ is 
the NSI-induced effective rotation angle in matter. Both 
the effective mixing angle and the correction factor de- 
pend on the neutrino energy as well as the standard oscil- 
lation and NSI parameters. These parameters are shown 
as a function of energy for several values of e and e' in 

Fig.m 

B. Expected Phenomena 

Because of the energy dependence of the standard os- 
cillation term in Eq. ^ , the relative dominance of NSI in 
the hybrid model is expected to depend on the neutrino 
energy. The effects of NSI on neutrino oscillation can 
be divided into three energy ranges: (1) E,^ <lGeV, (2) 
1< <30GeV, (3) E^ >30GeV. 

(1) E^ <lGeV 

At these low energies, the eigenvalue of the vacuum 
term Am23/2E^{> 1 x 10~^^eV) is larger than that of the 
NSI matter term •\/2G'F-/V/e(~ 1 x 10~"^'^eV), assuming 
Amia = 2.1 x lO^^cV^, Nf = Na ^ 3Ne and e - 0(1). 
Thus the i/^ — > i/t- transition is mostly governed by the 
standard two-flavor oscillation and there is no significant 
contribution from NSI. Note that e ^ 0{1) is a conser- 
vative assumption according to the NuTeV limit le^rl < 

0.05 mi. 

(2) 1< E^ <30GeV 

This energy range corresponds to the FC Multi-GeV 
(visible energy greater than 1330 MeV), PC, and UPMU 
stopping samples. In this region the matter term com- 
petes with the vacuum term and thus the Vfj, — ^ Vr transi- 
tion is no longer dominated by standard oscillations but 
is modified by the matter term. The left panel of Fig. [T] 
shows that for nonzero NU the effective mixing angle 
decreases with increasing neutrino energy, thereby sup- 
pressing — >■ transitions. Similarly, the right panel 
shows that nonzero FCNC interactions affect the fre- 
quency of oscillations. Since the effective mass splitting, 
Am^^ = i?Am23, is larger than Am^g, the first oscilla- 
tion maximum is expected to occur at higher neutrino 
energies than for standard oscillations. Focusing on the 
zenith angle distributions presented in the left panel of 
Fig. [21 the magnitude of the //-like deficit in the upward- 



going direction is expected to become smaller due to e'. 
Further, the shape of the zenith angle distribution near 
and above the horizon is modified in the presence of e for 
the higher energy samples. 

(3) E^ >30GeV 

Above 30GeV, most atmospheric neurinos are in the 
UPMU through-going sample, which ranges in energy 
from tens of GeV to ^ 100 TeV.At these energies vacuum 
oscillations have less of an effect on the i/^ — >■ v-r tran- 
sition, while the transitions induced by the matter term 
become significant when the neutrino path length in the 
Earth is sufficiently long. As shown in the right panel 
of Fig. O FCNC interactions play a leading role in this 
energy range because the modified oscillation frequency 
is comparable in size to the oscillation frequency of these 
neutrinos in the Earth. Thus, — > Vr transitions driven 
by FCNC are expected to occur. In contrast, at these en- 
ergies NU interactions suppress the effective mixing angle 
and therefore are expected to produce a sub-leading ef- 
fect on the data. Note that at the current values of the 
atmospheric mixing parameters, the standard z^^ — ^ Vr 
transition is already increasingly supressed at these en- 
ergies. In summary, FCNC interactions induce a faster 
oscillation frequency and lead to t'^ — >■ Vr transitions at 
shorter path lengths in matter, while NU suppresses the 
transition even at these energies. 



C. Analysis Method 

We evaluate the agreement between the data and 
Monte Carlo oscillated according to the hybrid NSI model 
using a test. The SK-I and SK-II data are divided ac- 
cording to their reconstructed event types, momenta, and 
zenith angles into 400 and 350 bins, respectively. In SK-I 
(SK-II) there are 310 (280) bins for the FC samples, 60 
(40) for the PC samples, and 30 (30) for the upward- 
going muon samples. Data from the two run periods are 
treated separately due to differences in the detector re- 
sponse and in the effects on the atmospheric neutrino 
flux from solar modulations during the runs. In order to 
accurately treat bins with small statistics in this bin ning 
scheme, a likelihood based on Poisson probabilities [2g 
is used. The complete with 750 data bins and 90 
systematic uncertainties is defined as 



750 



90 



obs 



90 

E 



(4) 



where cj is a fitting parameter for the j-th systematic the i-ih bin due to a Icr change in j-th systematic error, 
error and /j is the fractional change of the event rate in 




FIG. 1: (Left) Effective mixing angle in matter. (Right) Effective mass squared difference in matter. In botfi panels, solid black 
curves indicate the case with (e (FCNC), e' (NU)) = (0.015, 0.05), dashed red curves with (0.015, 0.0), dotted green curves with 



(0.0, 0.015), dashed-dotted blue curves with (0.0, 0.05). 
are assumed. Matter density is defined as constant p 
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FIG. 2: Zenith angle distribution for PC through-going (left) and UPMU showering (right, included in UPMU through-going 
where energy loss of muon is caused by pair production, Bremsstrahlung, and photonuclear interactions) sub-samples. In 
the solid green line (e, £') = (l-0 x 10"^, -2.4 x 10"^), in the dashed blue line (1.0 x 10~^, -0.38), and in the dotted red line 
(3.2 x 10"^, -2.4 x 10~^). In all lines the standard oscillation parameters are 023 = 45° and Amis = 2.1 x lO^^eV^. 



Using this equation, a value of is evaluated at each 
point in a four-dimensional parameter spaee defined by 
sin'^ 26*23, ^"i23, logioE, and logioe', where e and e' range 
from 1.0 x 10~^ to 3.2 x IQ-^ and from 1.0 x IQ-^ to 0.42, 
respectively. The fit is performed on a 51x51x51x51 
grid in this space. In Eq.(|4]), the 90 Cj parameters are 
varied to minimize the value of for each choice of os- 
cillation parameters. The point in parameter space re- 
turning the smallest value of is defined as the best fit. 
Since the effects of e are symmetric between negative and 
positive values, we only consider positive e. 



D. Two-Flavor Non-Standard Interaction Analysis 
Result 

The result of a scan on this parameter space gives a 
best fit at 

sin^ 26*23 = 1-00, Am^g 2.2 x 10~^eV^ 
e = 1.0xlO-'\ e' = -2.7xl0"2 

xL„ = 838.9 / 746.0 d.o.f. (5) 

The best fit value assuming only standard 2-fiavor os- 
cillations with the same samples and binning is Xmm = 
834.3/748 d.o.f. at sin^ 26* = 1.00 and Am^ = 2.1 x 
10~'^eV^. The sfightly larger Xmm value from the two- 
flavor hybrid model compared to the standard analysis 
is caused by the different numerical method used in the 
NSl analysis: the chosen ranges of e and e' do not al- 
low the two-ffavor hybrid model to reduce completely to 



standard oscillations. 

Fig. [3] shows the allowed neutrino oscillation parame- 
ter regions, sin^ 2^23 and Am^^, from this analysis and 
from the standard oscillation analysis. The three con- 
tours correspond to the 68%, 90% and 99%C.L., defined 
by = xliin +2.30, 4.61, and 9.21, respectively There 
is no inconsistency between the two sets of allowed re- 
gions. Since the difference in the minimum values is 
not large, no significant contribution to standard two- 
flavor oscillations from NSI effects is found in this anal- 
ysis. 
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FIG. 3: Allowed oscillation parameter regions derived by the 
2-fiavor hybrid model analysis (solid curves), where undis- 
played parameters e and e' are integrated out. For reference, 
the result of standard 2-flavor oscillation is added (dashed 
curves) . 
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The allowed regions of the NSI parameters are shown 
in Fig. m where there undisplayed parameters sin^ 26*23 
and Am23 have been minimized over. At 90% C.L. the 
obtained limits on the NSI parameters in the i'^ — 
sector arc 

|£| < 1.1 X IQ-^ and - 4.9 x 10"^ < e' < 4.9 x 10^^(6) 

These limits can be compared with a phenomenologi- 
cal study using the SK-I (79 kton yr) and MACRO atmo- 
spheric neutrino data [1^1 , in which the limits on FCNC 
and NU at 90% C.L. are 2.0 x 10"^ < e < 1.3 x lO'^ 
and -4.7 x IQ-^ < e' < 4.2 x lO'^. 



E. Discussion 

As described in Sec. IIVBI NSI are expected to affect 
multi-GeV muon neutrinos , and are consequently con- 
strained by the Super-K ^-like samples. Fig. [5] shows 
the allowed region spanned by the NSI parameters for 
three sub-samples of the data. Since the minimum 
is located in negative e' space, the negative e' plane is 



FIG. 4: Allowed NSI parameter regions in the 2-flavor hy- 
brid model. The horizontal axis shows e (FCNC) and the 
vertical axis shows e' (NU). Undisplayed parameters sin^ 2623 
and Am23 are integrated out. The star represents the best-fit 
point for the NSI parameters. 



presented. The strongest constraint on s comes from the 
UPMU through-going sample (solid curve), while e' is 
most tightly constrained by the PC and UPMU stopping 
(dashed curve) samples. The constraint from the sub- 
GeV samples is too weak to be visible in this figure. 

These constraints on e can be understood as follows. 
As shown in the right panel of Fig. [TJ FCNC interactions 
(dashed red curve) increase the effective neutrino oscilla- 
tion frequency thereby shortening the oscillation length 
for neutrinos at a fixed energy. Above 30 GeV, where 
the atmospheric sample is dominated by UPMU events, 
the oscillation length in the absence of NSI is already 
greater than the diameter of the Earth. The addition of 
FCNC effects on the other hand can shorten the oscilla- 
tion length enough to induce oscillations in the steepest 
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FIG. 5: Allowed NSI parameter regions for various event sub- 
samples. The solid curve indicates the allowed region from the 
UPMU through-going sample, the dashed curve is that from 
the PC and UPMU stopping samples, and the dotted curve is 
from the PC Single-ring Multi-GeV and Multi-ring samples. 
Each contour corresponds to = + 2-30 (68%C.L.). 



upward-going bins of the UPMU sample. This is an ap- 
parent contradiction with the data and results in a tight 
constraint on e. Below ^ 30 GeV, the oscillation length 
is sufficiently short in both standard and NSI cases to 
cause oscillations in the FC and PC samples. For this 
reason it is more difficult for these samples to discrimi- 
nate between standard- and FCNC-induced oscillations. 

As for the most stringent limit on NU derived from the 
PC and UPMU stopping samples, since a large e' gives a 
small magnitude of the //-like deficit in the upward-going 
direction and a distortion of the shape of the zenith angle 
distribution near and above the horizon (see dashed blue 
line in the left panel of Fig. ^ , a tight constraint on NU 
can occur. Further limits on NU by other sub-samples 
can be understood by the discussion in Sec. IIVBI 

V. ANALYSIS WITH A THREE-FLAVOR 
HYBRID MODEL 

We now consider a 3-flavor hybrid model in which 
NSI in the v^. — Vr sector coexist with standard 2-flavor 
O- Vt oscillations and all other NSI are zero. By 
introducing couplings between and this model al- 
lows flavor transitions of all types. That is, an overall 
j/^ — >■ j/g transition becomes possible due to the the £er 
induced Vt — > Vf, conversion working in conjunction with 
the standard i/^ -;-> v-^ oscillation: 

Vpi > Vr > Ve- (7) 

Note that since the other FCNC cpsilons have been set 
to zero, there is no direct transition between and u^. 
This analysis therefore aims to constrain possible NSI 



effects by examining changes in the and j/^ fluxes on 
top of standard oscillation-like effects. 



A. Formalism 

An evolution matrix from time to t can be obtained 
by diagonalizing the Hamiltonian Eq. ([1]) in terms of the 
leptonic mixing matrix in matter U' and the effective 
eigenvalues H ~ diag(i?i, i?2, In the case of constant 
matter density, the evolution matrix in natural units is 

3 

i=l 

(8) 

Thus the neutrino oscillation probability under the effect 
of NSI can be expressed as 

= |^/3a(t,to)r (9) 

However, in order to account for the varying matter den- 
sity in the Earth (from 2.5 to 13 g/cm'^), we divide the 
neutrino path into several constant density steps and cal- 
culate the evolution matrix in each. The oscillation prob- 
ability is obtained from the eigenvalues of the product of 
these matrices. 



B. Expected Phenomena 

As in the 2-flavor hybrid model, the effects of NSI are 
expected to vary with the neutrino energy. We consider 
three energy ranges: (1) E^, <lGeV, (2) 1< <15GeV, 
and (3) E^ >15GeV. 

(1) E^ <lGeV 

Since the standard oscillation eigenvalue lS.m? /2Ei, is 
much greater than the matter potential ^j^GrNd, — > 
Vr transitions induced by these oscillations are expected 
to be dominant and the effects from NSI can be ignored. 

(2) 1< E^ <15GcV 

In this energy range, the matter term has a sizable 
effect on the Vfj, — > Vr transition as it did in the 2-flavor 
hybrid model. Moreover, since the current limits on the 
NSI parameters governing Vf, — >■ Vr transitions are poor: 
Ser ~ 0(10"^) [201, large effects from NSI are possible. 

The modification of the Ve flux by NSI can be param- 
eterized by two oscillation probabilities: P(;/e — Ve) and 
P(i-'p Vg) in Eq. Normalizing by the flux($e), 
and disregarding any contamination, the resulting 
e-like distributions can be approximately expressed as 
P{Vf. -> Ve) +rP(vLi — » Ve)^ where r = $^/$e is the neu- 
trino flavor ratio [27|. This ratio grows with increasing 
neutrino energy and though it is highly dependent on the 
neutrino zenith angle, near the horizon it is ^ 2 up to 
around 10 GeV. Thus, if P{ve Ve)+rP{v^ ^ Ve) < 1, 
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as is the case when Ser and Stt have comparable val- 
ues, then since P{ve v,,) is suppressed by £er and 
P{^^l ^ i^e) by Srr, the number of e-like events at the 
horizon can be expected to decrease. In contrast, the 
number of e-like events in the upward-going bins would 
effectively increase because large values of Ser produce 
additional v^. from Vr created during standard oscilla- 
tions. This effect is enhanced by the larger flavor ratio in 
these regions and enables the multi-GeV e-like samples 
to help constrain e^T ■ These effects are shown in the left 
panel of Fig. [51 

Effects driven by Srr on the other hand may be de- 
scribed in terms of the limiting case where the other NSI 
are set to zero. In this case the problem reduces to the 
2-flavor hybrid model with e = and e' = Ett- Therefore 
a constraint on e' similar to that from the 2-flavor case 
can be expected. 

(3) >15GeV 

Above a few tens of GeV the I'e flux decreases as their 
parent muons increasingly reach the ground before de- 
caying, making it possible to neglect the Ve contribution 
to the h'fi flux from oscillations induced by Ser- Con- 
versely, changes in that flux can be clearly recognized as 
NSI-driven i/^- — > ly^ conversion at these energies. As in 
the 2-flavor hybrid be expected to suppress 

the effective mixing angle and increase the effective mass 
splitting in the z^^ — s- i^r sector. Both of these effects can 
be seen in the right panel of Fig. |6l 

As pointed out in Ref. j29| the atmospheric neutrino 
sample cannot effectively constrain See- This is due to the 
fact that when Ser is zero, the eigenstates of the Hamil- 
tonian in matter are identical to the vacuum eigenstates. 
In this case, the matter eigenvalues are no longer depen- 
dent upon Eee and the problem reduces to two-flavor NSI 
mixing with i/f^ <-> h'r transitions in matter modified by 
only Stt- Despite the lack of sensitivity to See, at fixed 
non-zero values it dictates a parabolic relationship among 
the NSI parameters, which in the limit where one of the 
matter eigenvalues is small takes the form 



This is a general feature which will present itself in the 
shape of our allowed regions below. 

C. Analysis Method 

The analysis procedure for the 3-flavor hybrid model 
follows that used for the 2-flavor hybrid model. A 
value of is evaluated at each grid point in the three- 
dimensional parameter space of See, Eer and Ett, where 
51 points are chosen for each. Since Est enters the oscilla- 
tion equations as |eeT-| [29j when 813 and solar oscillations 
(Am^2 and 612) are not considered, we only test positive 
values of the parameter. 





Best-fit Est 


Best-fit Ett 


Minimum x'^ 


-0.50 


0.016 


-0.016 


831.1 


-0.25 


0.016 


0.024 


829.9 


0.00 


0.024 


0.016 


830.9 


0.25 


0.000 


-0.016 


831.4 


0.50 


0.000 


-0.016 


831.4 



TABLE I: The allowed NSI parameters at the 90% C.L. as a 

function of See- 



We first set the standard 2-flavor parameters to 
(sin^6'23, Ato2)=(0.5, 2.1xlO"3eV2) as motivated by the 
results from the 2-flavor hybrid model and other Super-K 
analyses [l^. Here the standard oscillation parameters 
are taken as fixed values. This assumption will be ver- 
ified later by comparing the allowed NSI parameter re- 
gions derived using slightly different standard oscillation 
parameters. Modifications to the fitting results when ^13 
is non-zero are presented in Sec. I VI I 



D. Results of the 3-Flavor NSI analysis 

The allowed regions for the 3-flavor NSI parameters 
are shown in Fig. [T] Allowed regions for five fixed values 
of See (-0.50, -0.25, 0.00, 0.25, and 0.50) are presented. 
The three contours correspond to the 68%, 90% and 99% 
C.L. regions as defined by = Xmm +2-30, 4.61, and 
9.21, respectively. The best-fit values for each See are 
summarized in Tbl. HI 

In these figures, except for the projection for See = 
—0.25, parabola-like regions in the Eer - £tt are shown. 
In particular, the allowed region extends to negative Ett 
values when Eee = —0.5, while they extend to the positive 
values above E^e = 0. This feature is consistent with a 
transition of the matter eigenvalue hierarchy discussed in 
Ref. 

Next, we evaluate the influence of using flxed standard 
oscillation parameters. Fig. [5] shows the allowed NSI pa- 
rameter regions using four slightly different sets of stan- 
dard oscillation parameters. The oscillation parameters 
are taken from the 90% C.L. allowed region from the 
standard SK 2-flavor analysis. Although slight changes 
appear when the NSI fits are repeated using these param- 
eters, their minimum values are larger than that of the 
original fit by less than 1.5 units. Therefore no significant 
change in the allowed parameter regions is expected even 
if the standard oscillation parameters are allowed to vary 
during the NSI fit. Finally note that the current allowed 
regions of sin^ 623 and Am^ are constrained more tightly 
by analyses using the full SK data set, which is larger 
than the sample used here. 
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FC Multi-GeV 1-rinq e-like 




0.8 0.6 0.4 0.2 0.2 0.4 0.6 0.8 1 

cos 6 



PC throughgoing 




0.8 0.6 0.4 0.2 0.2 0.4 0.6 0.8 1 

cos 6 



FIG. 6: Zenith angle distribution for typical e-like (left, FC Multi-GeV 1-ring e-like) and /x-like (right, PC throughgoing) 
samples. The solid green line is the MC prediction at See ~ 0.0, Ser = 0.2, and Ett = 0.2. The dashed red line is that for 
standard neutrino oscillations. In all lines the standard oscillation parameters are 623 = 45° and Am23 — 2.1 x 10"'' eV^. 



<o^ ; SK-i + SK- 

0.15 - 




99% C.L. 
90% C.L. 
68% C.L. 



(0^ ; SK-I + SK- 

0.15 - 



99% C.L. 
90% C.L. 
68% C.L. 



hO.2 
0.15 



; sK-i + sK-ii 






99% C.L 




90% C.L. 


- 


68% C.L. 



(a)ee 



-0.50 



ih)e 



-0.25 



(c)£e 



: 0.00 



CO 

0.15 



;sK-i + SK-ii 






99% C.L. 




90% C.L. 


- 


68% C.L. 



0.15 



;sK-i + SK-ii 






99% C.L. 




90% C.L. 


- 


68% C.L. 



(d)e 



: 0.25 



(e)£e 



: 0.50 



FIG. 7: Allowed NSI parameter regions at fixed Eee at 68%, 90% and 99% C.L., where contours are drawn at — Xmin = 
2.30,4.61, and 9.21, respectively. 



E. Discussion 



Fig. [5] shows the allowed parameter regions for three 
subsets of the atmospheric sample. The high energy v^- 
rich UPMU throughgoing sample (indicated by the solid 
curve) and PC-I-UPMU stopping samples constrain Err 
as expected. However, these samples provide no signifi- 
cant contribution to the constraint on Ser because they 
lack an e-like component. On the other hand, the FC 



multi-GeV samples (shown by the dotted curve) include 
several e-like sub-samples which better constrain Ser- 



VI. ANALYSIS WITH A THREE-FLAVOR 
HYBRID MODEL WITH NONZERO On 

There are several scenarios in which sub-dominant ef- 
fects from standard neutrino oscillations may affect the 
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"0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 p 0.4 



(a)ee 



-0.50 



(b)£e 



: 0.00 



(c)£e 



: 0.50 



FIG. 8: Allowed NSI parameter regions in the Est vs. Ett plane at the 90% C.L. for Eec ~ -0.5, 0.0, and 0.5 using four sets of 
standard oscillation parameters. The solid curve corresponds to (sin^ 623, Am^) — {0.5, 1.7x 10~'^eV'^), the dashed curve to (0.5, 
2.7xlO"^eV^), the dotted curve to (0.39, 2.1x 10"^eV^), and the dashed-dotted curve to (0.61, 2.1xlO"^eV^). 




£0.4 



FIG. 9: Allowed NSI parameter regions at 90% C.L. derived 
for three subsamples at Eee = —0.25. The solid curve indicates 
the allowed region given by the UPMU throughgoing sample, 
the dashed curve shows the PC and UPMU stopping sample, 
and the dotted curve is the constraint from the FC multi-GeV 
samples. 



allowed NSI parameters. In this section we consider how 
our limits change when effects from nonzero ^13 are in- 
cluded in the analysis. 

Present data suggest that 6*13 is small relative to the 
other mixing angles, 623 and 9i2. The Chooz experiment 
has placed the most stringent limit on the parameter, 
indicating that sin^6'i3 < 0.04 [s^. Due to its small size 
its effects were ignored in the main analysis. However, 9i3 
is expected to induce t'^ — >■ Ve transitions at multi-GeV 
energies possibly producing an excess of that could be 
misinterpreted as the effect of Ser- Here the analysis of 
section|V]is repeated for the normal hierarchy, ATO23 > 0, 
and inverted hierarchy, Am^^ < 0, with 613 fixed at the 
Chooz limit. Including ^13 breaks the positive-negative 



symmetry of Eer so the fit is expanded to cover both 
regions. The results of the fit are again presented as a 
function of See hi Fig. [TOl 

Focusing on the normal hierarchy case, the addition of 
9i3 to the fit tends to improve the constraint on Ser ■ The 
— > oscillation probability in constant density mat- 
ter is proportional to sin^0i3, where Q13 is the effective 
613 mixing angle given by. 



61 



tan 20 ■ 



a sin 29 



13 



At- 



''31 



a cos 29 



13 



(11) 



In this equation a is the product of the MSW matter po- 
tential and the neutrino energy, ±2^/2GFNeEy, where 
the sign is positive (negative) for neutrinos (antineutri- 
nos) . The structure of the denominator can create a res- 
onant enhancement of the oscillation probability depend- 
ing upon the energy of the neutrino and the density of 
matter it traverses. If 9i3 is nonzero, the upward-going Vf. 
flux is expected to increase in the 2-10 GeV range, coin- 
ciding with the region where Ser can enhance the v,, flux. 
Accordingly, part of the Ve flux induced by Ser is now 
occupied by events from 9i3 transitions which therefore 
results in a tighter constraint on the parameter. How- 
ever, the resonance behavior of Eqn. (|lip is contingent 
upon the signs of the mass hierarchy and MSW matter 
potential. This, coupled with unequal neutrino and an- 
tineutrino fluxes in the atmospheric data, results in the 
asymmetric constraint on Ser seen in Fig. 1101 



VII. CONCLUSIONS 

We have studied non-standard neutrino interactions in 
the context of atmospheric neutrinos propagating in the 
Earth. Two analyses were presented considering possi- 
ble effects from both flavor changing neutral current and 
lepton universality violating interactions. Analysis of the 
SK-1 and SK-11 atmospheric neutrino data shows no evi- 
dence of NSI and provides the following constraints. For 
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(c)£e 



: 0.00 



(0 
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-0.05 
-0.1 
-0.15 

n. 



SK-I -I- SK-I 



0.15 



-0.15 



SK-I -I- SK-I 



(d)ee 



: 0.25 



(e)ee 



: 0.50 



FIG. 10: Allowed NSI parameters regions in the Ebt vs. £tt plane at the 90% C.L. when sin'^^is = 0.04. The solid curve and 
dashed curve indicate the normal inverted hierarchy fits, respectively. 



NSI in the v^ — Vt sector, the 2-flavor hybrid model allows 
contributions from NSI in the form 

|e^r| < 1-1 X 10"^ and 
-4.9 X 10-2 <£rr - £mm < 4-9 X 10^2 , (12) 

at the 90% C.L., where e and e' are replaced with e^r and 
Ett — £(1(1 1 respectively. In the 3-flavor hybrid model, the 
allowed regions are presented for different values of e^e 
since the atmospheric data have no ability to constrain 
this parameter. 
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